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The^paper  gives  the  general  theory  of  optimal  confidence  '  o'  *  « 

\*  i  y 

limits  for  systems  reliability  introduced  by  Buehler  (1957).  V  f J 
This  is  specialized  to  series  systems.  It  is  noted  that  some 
results  previously  given  are  false.  In  particular,  counter¬ 
examples  for  results  of  Sudakov  (1974),  Winterbottom  (1974) 
and  Harris  and  Soms  (1980,  1981)  are  given.  Numerical  examples 
are  provided,  which  suggest  that  despite  the  theoretical  problems 


of  the  results,  they  are  nevertheless  valid  for  significance 


levels  likely  to  be  used  in  practice 
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SIGNIFICANCE  AND  EXPLANATION 


Systems  with  independent  components  arise  naturally  in 
engineering  practice.  Therefore  it  is  of  importance  to  efficiently 
utilize  data  obtained  on  individual  components  to  obtain  an  assess¬ 
ment  of  the  reliability  of  the  system. 

This  paper  presents  a  unified  theory  for  doing  so  and  points 
out  errors  in  previous  results  on  series  systems. 


The  responsibility  for  the  wording  and  views  expressed  in  this 
descriptive  summary  lies  with  MRC,  and  not  with  the  authors  of  this 
report . 
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Bernard  Harris  and  Andrew  P.  Soms 

1.  Introduction  and  Summary 

A  problem  of  substantial  importance  to  practitioners  in 
reliability  is  the  statistical  estimation  of  the  reliability  of  a 
system  of  stochastically  independent  components  using  experimental 
data  collected  on  the  individual  components.  In  the  situations 
discussed  in  this  paper,  the  component  data  consist  of  a  sequence 
of  Bernoulli  trials.  Thus,  for  component  i,  i*l,2,...,k,  the 
data  is  the  pair  (n^,Y^),  where  n^  is  the  number  of  trials  and  Y^ 
is  the  number  of  observations  for  which  the  component  functions. 

Yj , Yj , . . . , Y^  are  assumed  to  be  mutually  independent  random 
variables. 

This  problem  was  treated  in  Sudakov  (1974),  Winterbottom 
(1974),  and  Harris  and  Soms  (1980,1981);  one  purpose  of  the  pre¬ 
sent  paper  is  to  exhibit  counterexamples  to  theorems  in  the  above 
papers . 

In  Section  2  we  discuss  the  general  theory  of  optimal 
confidence  limits  for  system  reliability  so  that  the  notation  and 
definitions  to  be  employed  in  the  balance  of  the  paper  have  been 
prescribed . 

In  Section  3  the  counterexamples  previously  mentioned  are 
exhibited  and  the  specific  errors  in  the  proofs  of  the  theorems 
are  indicated. 
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Section  4  presents  the  proof  of  a  special  case  of  the  key 
test  theorem  (Winterbottom  (1974)),  the  general  form  of  which  was 
invalidated  by  a  counterexample  in  Section  3. 

The  consequences  for  reliability  applications  are  discussed 
in  Section  5. 


2.  Buehler's  Method  for  Optimal  Lower  Confidence 
Limits  for  System  Reliability 


We  now  introduce  the  notation,  definitions,  and  assumptions 
that  will  be  us~d  throughout  the  balance  of  this  paper, 

1.  Let  p^,  i«l,2,...,k  denote  the  probability  that  the  ith 

component  functions.  The  components  will  be  assumed  to  be 
stochastically  independent.  The  reliability  of  the  system  will 
be  denoted  by  h(p),  where  p  ■  (Pj ,P2 » • • • »Pk) »  0  <  <  1.  It  is 

assumed  that  h(0,0,...,0)  «  0,  h(l,l,...,l)  -  1,  and  that  h(p)  is 
non-decreasing  in  each  p^,  i«l,2,...,k.  Further,  h(p)  is 
continuous  on  {p|0  <  p^  <  1),  which  follows  readily  from  the 
assumption  of  independence.  These  properties  hold  for  coherent 
systems  (see  Barlow  and  Proschan  (1975)). 

2.  Let  S  -  {i|»i  -  0,  L,...,ni#  , 2 , . . . , k }  be  the  failure  set. 

g(x)  is  said  to  be  an  ordering  function  if  for  x^  <_ 

x2  <_  *2,...,xk  <  z  k ,  x,  z  e  S,  g  lx)  >_  gtz).  (It  is  often  convenient 

to  normalize  g(x)  by  letting  g(0)  ■  1  and  g(n)  ■  0.  With  such 

a  normalization,  glx)  is  often  selected  to  be  a  point  estimator 
of  h(p).) 

3.  Let  R  ■  {rj ,r2, . . . ,rs ,  s  >  2)  be  the  range  set  of  g(x),  With 

no  loss  of  generality  we  order  R  so  that  Tj  >  r2  >  . , .  >  rs< 


-2- 


■  ^  -  •  _  •  .  *  _■  _  *  _  a  _  r  .  -'a*  *  *  * 


« 


We  now  introduce  the  notion  of  Buehler  optimal  confidence 


limits.  Let  g(x)  ■  r j .  Then  define 


*St5)  '  inf 


{h(P)|Pf{i| 


8(i)  >  8(x) 


M- 


(2.5) 


Equivalently,  by  (2.2),  we  can  also  write 


ag  (S) 


inf|h(p)|P~|xe  U  aJ  >  aj-  . 


(2.6) 


We  now  establish  the  following  theorem. 


Theorem  2.1.  Let  assumptions  1-5  be  satisfied.  Then,  for  xeS, 

a^^jj  is  a  l-o  lower  confidence  limit  for  h(p).  If  bg^  is  »ny 

other  1-a  lower  confidence  limit  for  h(p)  with  b  >  b_  >  ...  >  b  , 

1  ’  2  ‘  ’  ) 

then  <  ag(x)  for  a11  *eS* 

Proof.  Fix  p  and  let  m(p)  be  the  smallest  integer  such  that 


“(P)  ) 


>  o  . 


Then 


P-(x  e  U  A, 1  >  l-o  . 

Pv  i«m(p)  * 

\  ■  {5|pp{;  c  j,  *i}  i  “}  • 


Then  D  is  a  1-a  confidence  set  for  p ,  since 
g(X) 


*D«C*)}  ’  Pi{<"'x>  i  r»(J)}  i  l-a  • 


•*.  •*.  «■-  **.  .  •  •  -  ••  •  *  •.*  •.*  %•  *,•  ■;  •.  •:  ;  \v». /•  /-V-  /•.** v 


By  assumption  1,  h(p)  is  continuous  and  the  set  of  parameter 
points  satisfying  (2.5)  is  compact;  therefore  the  infimum  in  (2.5) 
and  (2.6)  is  attained. 


Assume  that  there  is  an  integer  j,  1  <  j  <  s-1,  such  that 
br  >  ar  .  Then  there  exists  a  pQ  such  that 

br^  >  a^  -  inf|h(p)|P^|x  e  aJ  >  o}  -  h(pQ)  .  (2.7) 


In  addition,  there  exists  a  p ^  such  that 


Pr  jxe  u  A. \ 

Pll  i.i  ii 


>  «  »  h(P.j)  K  br 


(2.8) 


Since  b  >  b  >  ...  >  b  ,  from  (2.7)  we  have 
1  ‘  2  ‘  s 


b(p.)  <  b  ,  t  »  1 , 2, . . .  ,  j  . 
1  Tt 


Therefore 


(2.9) 


'Pxf* £  j,  *i}  ■  'p,^’  ‘  b«m}*  <2-10) 


which  is  a  contradiction.  Consequently,  there  is  no  integer  j. 


i  <  j  <  s-1,  for  which  b„  >  a,  , 

•  •  rj  ri 

From  (2.6),  it  follows  that  a  *0  and  br  is  also 

s  s 


necessarily  zero.  Note  further  that  in  (2.7)  it 
that  the  infimum  is  attained  at  a  point  for  which 
To  see  this  consider  the  following  example. 


is  possible 

ppCx  c  L  ai} 


> 


Let  k  »  2,  nj  »  5,  n2  ■  10,000,  Xj  -  0,  x2  *  5, 
g(x)  ■  nj+nj-Xj-Xj,  h(p)  ■  PjP2.  11  *s  seen  that  the 

hypotheses  of  Theorem  2.1  are  satisfied.  Thus,  for  the  data 


given,  g(x)  «  10,000  -  r The  set  U  A^  consists  of  all  points 
(xlfx2)  for  which  Xj+x2  <  5,  that  is,  Aj  -  |(0,0)j, 

*2  "  |ci»0) , (0, 1)J,  and  *o  on.  Consequently, 

includes  the  parameter  points  (0,  p2a)  where  p2a  satisfies 

Pp2o{X2“0}  ^  ®*  since  Pjj{xi  <  s|  -  1  when  Pj  -  0.  Thus  inf 
h(?)  ■  0  for  all  0  <  a  <  1. 

We  note  that  the  nonotonicity  of  h(p)  is  not  utilized  in 
the  proof,  which  is  valid  whenever  h(p)  is  continuous. 

It  is  easy  to  see  that  is  ■onotone,  i.e.,  ay  >  ar  > 

This  follows  from  (2.7)  upon  noting  that  as  j 
s 

increases,  the  set  of  p  satisfying  (2.7)  increases  and  the 
infinum  is  taken  over  a  larger  set. 

k 

Corollary.  For  a  series  systea  h(p)  »  Up..  Then  if 
“  k  k  i-1  1 

g(x)  -  H  (n.-x.)/n.  ■  II  y./n4,  the  hypotheses  of  Theoren  2.1 
i-1  i-1  1  1 

are  satisfied  and  the  conclusion  follows. 

k 

Note  that  g(x)  ■  n  (n.-x  )/n.  is  the  maximun  likelihood 

i-1  11 

estimator  as  well  as  the  miniaun  variance  unbiased  estimator  of 
k 

II  p.  and  is  therefore  a  reasonable  choice  of  an  ordering  function 
i-1 


We  now  establish  the  following  theorem. 


Theorem  2.2.  Let  g(x)  -  r^  and  let 


f  (x;a)  -  sup 
h(p) 


,t  ;  *,}.  o 


<  a  <  1  . 


(2.10) 


s“  •  .  *  .  *  *  .‘s'  *  V  '  \r\*t  •  .‘X'T  %  *  .  O”  S*V*  •  O’ /  v‘. 

>  .  •  .  *,».•.*  ,  ».•  ,  «  ,*  .  1  •  ,  •  ,  *  ,•  ,  *  .  ‘j,  ^  *  *  .  »  ,  '>>>  i  %  ^  ^  v  %  '  *  -  ■ 

*  - - 


I  .«*riL^l>i4CUK.  * 


Then 


sup  f*(x;a)  ■  1 
0<*<1 


and  f  (x;a)  is  non-decreasing  in  a. 

Proof.  Since  h(p)  is  continuous  and  h(l)  »  1, 


lie  sup  P=fg(X)  >  r,j  »  1  . 
.-►I  hfn'i-a 


a-*-l  h(p)»a  **  ' 

Now  choose  a  and  b  such  that  0<a<b<l, 


*}{zm  >  r, 


J  •  f*(x;j 


“/V 


pp  {*CXJ  >  rj  -  f#(x;b)  . 


Let  Ift  be  the  set  of  indices  i  such  that  p^a  <  1.  Then  it  is 
possible  to  replace  pia  by  P^b,  ie  I#,  where  pia  <  p£b  <  1,  so 
that  h(p^)  «  b,  where  p^b  ■  Pia»  ie  I*.  This  follows  since 
h(l)  »  1  >  a  and  h(p)  is  continuous.  The  conclusion  follows  fron 
the  monotone  likelihood  ratio  property  of  the  binomial  distribution. 

Note  again  that  only  the  continuity  of  h (p)  was  used  in 
the  proof  of  Theorem  2.2. 

For  the  case  of  series  systems,  it  is  possible  to  strengthen 
Theorea  2.2  and  to  exhibit  the  above  construction.  This  is  done 
below. 

k 

Corollary.  Let  g(x)  »  r  .  If  h(p)  ■  JI  p.  ,  then  inf  f  (x;a)  »0 
#  *  i»l  0<a<l 

and  f  (x;a)  is  strictly  increasing  in  a  whenever  all  u^  <  n^ 

(see  (2.3)  for  the  definition  of  u^),  j«l,2,...k. 


,  V 

«  *•  1 

,  -v<, 

V*  «*• 


vy-  :*• 


-J 

•A-  *>.H 
y.v. 

-.t  -j 


•**  w*  *«*  .*  j*  _<•  \  ’  »  '  •  L  »  '  U  ■  .  *  •  •_#»»«,  ••*»'*.  •  •**«**■-•»  *  .  *  _  •  c  •  .  »  *  •  -  l_N 


Proof.  From  the  hypotheses 


"p{«m  >  *,}  <.  ‘V1  •  *■»•* . 1  • 

k 

and  since  n  p.  -*•  0  implies  at  least  one  p,  0,  this  gives 
i-1  1  1 


inf  f*(x;a)  *  0  . 
0<a<l 


To  show  that 
0<a<b<l  and  let  p 
Similarly,  let  p^ 

’•  *  {il,i2”  ' 
pai^a^1/r  <  1  Cn 


♦  -W 

f  (x;a)  is  strictly  increasing  in  a,  consider 

a  =  (pal»*“*pak5  satisfy  **(*;»)  *  Pp  |g(X)  > 
satisfy  f*(x;b)  *  P~  |g(X)  >  r.i.  Let 

be  any  non-empty  set  of  indices  such  that 
on-empty  because  otherwise  multiplying  the 


r 


(J 

components  would  give  b  >  1,  a  contradiction)  and  let  Ia  be  the 
remaining  indices.  Then 


1/r 


C  n  p  (— )  )  n  P 

jel  aIj  a  jei®  aij 


b  .  (2.11) 

'a  *  '  "a 

From  the  monotone'  likelihood  ratio  property  of  the  binomial 
distribution , 

PPa{gC*J  *  rj}  <  Pp*{gC*)  -  rj}  * 
where  the  components  of  p*  are  given  by  (2.11).  This  gives 


f*(x;a)  <  f*(£;b)  , 


which  is  the  desired  conclusion 


Note  that  if  at  least  one  u.  *  n  ,  it  follows  immediately 

3  3  k 

from  (2.5)  that  ag(£)  "  For  9(*)  =*  H  (n^-x^)/n^  the 

condition  <  n^  is  equivalent  to  x^  <  n  ,  j«l,2,...,k. 


/  J 


We  now  establish  a  result  which  will  prove  useful  in  some 
of  the  subsequent  material. 

Theorem  2.3.  If  f*(x;a)  ■  a,  0<a<l>  has  at  least  one  solution  in 
a,  then 

*g(x)  *  -  a}.  (2.12) 

If  f*(x;a)  >  a  for  all  a,  then  a...  ■  0. 

2  (x) 

Proof .  Let 

c  »  inf ja|f*(x;a)  >  aj  .  (2.13) 

The  infimum  in  (2.13)  is  attained.  Thus,  there  exists  a  pQ  such 
«.  *  ~ 

that  c  *  h(pQ).  If  f  (x;a)  >  a  for  all  a,  let  pt  0,  i*l,2,.... 
Then  h(p)  -*•  0,  since  h(0)  ■  0  and  h(p)  is  continuous,  and 

m«(i)  •  °- 

^  •** 

Now  assume  there  is  at  least  one  a  with  f  (x;a)  =  a.  Then 
f*(x;agj~j)  >  o  and  therefore  c  <  a^^-j.  If  c  <  ag ^J) ’  then 
c  ■  h(pQ)  and  f*(x;c)  =  a,  which  is  a  contradiction. 

Again,  only  the  continuity  of  h(p)  was  used  in  the  proof 
of  Theorem  2.3.  Under  the  hypotheses  of  the  Corollary  to 
Theorem  2.2,  for  a  series  system,  ag(jj)  i®  t^'e  s°iuti°n  in  a  °£ 

f*(x;a)  ■  a  .  (2.14) 

The  general  theory  described  in  this  section  applies  as  well 
to  what  is  known  as  systems  with  repeated  components  (see,  e.g., 
Harris  and  Soms  (1373)).  For  such  systems,  there  are  1  <  m  <  k 
unknown  parameters  P i  >?2  ’  *  *  *  'Pm*  s*nce  th®  "repeated  components'* 


•re  assumed  to  have  identical  failure  probabilities.  This 
assumption  permits  the  experimenter  to  regard  the  data  as 
^ni**i^’  i«l,2,...,m,  and  employ  the  previous  results. 

For  example,  if  a  series  system  of  k  components  has  otj  of 
one  type,  c*2  of  a  second,  ...,  of  an  m**1  type,  then 

Oj  a2  am  k 

MP)  »  p.  p,  . . .  p  ,  l  a,  -  k  . 

1  1  m  i-1  1 


3.  Counterexamples 


In  this  section  we  restrict  attention  to  series  systems  and 
employ  the  ordering  function 


g(x)  «  II  (n.-x.)/n.  , 
i®  1  11  1 


introduced  following  Theorem  2.1.  As  noted  previously,  in  this 

k 


case  the  reliability  function  h(p)  »  II  p.  .  With  this  special- 

i»l  1 


ization  we  have  for  (2.4) 


tl  * 


(3.1 


and  for  each  fixed  0  <  ij  <  t^,  0  <  i2  <  t2,  ...,  0  <  i^^  <  t^ 

j-1 


tj  «  (l*rm/C^n^  (nJl-iJl)/njl3) ,  2  <  j  <  k  , 


(3.2 


whenever  g (x)  ®  rffl,  1  <  m  <  s.  If  a  *  s,  then  rs  ■  0  and  aQ  *  0 


For  <  >  0,  X  >  0,  let 


ip(K,X)  -  r(^-Xy  |P  tK’1(l-t)X'1dt  ,  0  <  p  <  1  ,  (3.3: 


the  Incomplete  beta  function. 


*o  I 


It  is  well-known  that  if  t  is  an  integer,  t  <  n ,  we  have 


l  (?)  pn'V  ■  I  (n-t,fl)  .  (3.4) 

i-0  1  p 

Sudakov  (1974)  published  the  inequality 

p-{g(X)  >  rj  <  I  k  (ni-tj.t^l)  .  (3.5) 

n  P 

i-l  1 

This  inequality  and  generalisations  of  it  were  further 
studied  in  Harris  and  Sons  (1980,1981).  (3.5)  implies 

f*(x;a)  <  Ia(n1-t1,t1+l)  , 

hence  its  usefulness.  However,  as  we  now  establish,  (3.S)  is  not 
universally  valid,  as  was  claimed  in  Sudakov  (1974). 

Let  (Xj.Xj)  »  (Xj.O)  and  let  (n^nj)  ■  (n1,2n1).  Then 
g(x)  «  (n j -Xj)/nj  and  tj  ■  Xj .  Consider  P-|g(X)  >  rB|.  If 

■  (l,a),  0  <  a  <  1*  we  have 

Pp{««>  :  '»}  •  >  r.}  . 

since  P^X^ol  ■  1,  by  (2.1).  Consequently, 

V{*«>  i  >}.-*/.{**  -  n2C1-r»>} 

•  P.{x2  <  2»ltl-rn)}  . 

Since  rB  »  (n1-x]L)/n1, 

I  *.}  •  P.{*2  i  ixl)  ■ 


Thus  from  (3.4), 


Pp{g(X)  >  rn}  “  xa^2 (ni"xi^ *  2xi+1i  • 


The  Sudakov  inequality  implies  that 


Ia(2(nl-Xl)*  2xi+1^  -  Ia(nl-Xl»  xi+1) 


XaC2nlV  2nl(1"rm)  +  1)  -  Ia^nlrm*  nl(1-rm)+1)  *  (5‘6) 

Let  *12^t’n2*rm^  an<*  ^(tjn^.r^)  denote  th®  beta  density  functions 
corresponding  to  the  left  and  right  hand  side  of  (3.6),  respec¬ 
tively.  Then,  provided  nj^  >  1#  there  is  an  e  >  0  such  that 

h2(t;n2*rm)  <  hi(t;nl»rm)  0  <  t  <  e,  1-e  <  t  <  1  . 

This  implies  that  hjCt^j.r^)  and  h2(t;n2,rn)  intersect  in  at 
least  two  points.  If  t*  is  such  an  intersection,  setting 
h1(t;n1,rm)/h2(t;n2,rm)  *  1  gives 


n.  r  n, 

t  1  m(l-t)  1 


ni Cl-r  ) 


c(nx,rn)  >  0  . 


Thus,  for  1  <  m  <  s,  there  are  exactly  two  such  intersections. 


Therefore  there  is  a  z  such  that 

o 


lt  (nlV  ni^-rm)+1)  *  *z  Cn2rm*  n2(1'rm)  +  1)  * 


for  z  >  zQ, 


Iz^nlrm*  ni^-rm)+1)  <  Iz(n2rm*  n2(1“rm)+1) 


and  for  z  <  zQ, 


Iz(nlrm*,llC1*rm,  +  1)  31  Iz(n2rm'  n2(1“rm)  +  1) 


\  s 


Thus  for  z  >  z  .  (3.6)  is  violated.  (3.6)  was  used  as  a  lemma 
0 

by  Sudakov  (1974)  to  prove  the  inequality  (3.S).  This  lemma  was 
also  employed  in  Harris  and  Soms  (1980,  1981).  It  is  the  falsity 
of  this  lemma  which  invalidates  (3.5). 

Table  1  provides  some  illustrations  of  the  violation  of 
(3.5)  for  k  •  2  and  selected  values  of  (n1,n2)>  (x^,x2).  The 
smallest  value  of  P]P2  *0T  which  this  violation  occurs  is  also 
given  in  the  table,  where  it  is  denoted  by  a*.  In  addition, 

f*(x;a*)  is  tabulated.  Thus  for  o  <  f*(x;a*),  (3.5)  is  valid. 

The  calculations  were  made  by  means  of  a  FORTRAN  program.  Note 

the  for  (n1,n2>  -  (5,5)  and  (x^x^  -  (1,1)  the  inequality 

was  not  violated. 


Table  1.  The  Smallest  a,  a*,  and  f*(x;a*) 


(nlfn2) 

(x1.*2) 

a* 

£*(x;a*) 

(5,5) 

(1,1) 

1.0000 

1.0000 

(5,5) 

(3,3) 

.7454 

.9998 

(5.10) 

.  d.O) 

.8798 

.8909 

(5,15) 

(0.3) 

.8698 

.8791 

(5,30) 

(1.0) 

.8498 

.8467 

4  .  The  Theory  of 

Key  Test  Results 

If 

for  Uj  <  nj  <  ...  <  n^. 

(x j  ,  x2 ,  .  .  .  , Xj^)  *  (Xj 

,0, . . . ,0) , 

k  >  2,  then  x  is  called  a  key  test  result.  Winterbottom  (1974) 
asserted  that  subject  to  x^  <  f(k,n^),  where  f(k,n^)  is  the 


solution  in  f  of 


we  have  a  is  the  solution  in  a  of 

S 


Ia^nl"xl’  xi+1)  "  a  *  0  <  a  <  1  . 


(4.2) 


This  would  imply  the  inequality  (3.5),  which  we  have  disproved 
in  Section  3. 

As  we  subsequently  establish,  the  error  in  tfinterbottom's 
(1974)  result  is  a  consequence  of  falsely  concluding  that 
fChjn^)  depends  only  on  n^.  It  is  easy  to  be  led  to  this  conclu¬ 
sion  on  intuitive  grounds,  since  (n^-x^,  n^,...,  n^)  would  seen 
to  be  a  less  favorable  experinental  result  than  (n^-Xj ,n2, . . . ,nk) , 
whenever  for  at  least  one  index  i,  2  £  i  £  k.  We 

now  establish  a  modified  key  test  result  that  holds  for 
x^  <  f(k,n),  where  f(k,n)  is  the  smallest  solution  in  f  of 


k  -  „1/V 

n.  •  f  ■  1  •  k[(n. -f)  n  n.]1'*  . 
i-1  1  1  i-2  1 


(4.3) 


Theorem  4.1.  If  $  »2  f  •••  <  nk  and  *  *  (xx ,0, . . . ,0) ,  with 

Xj  <  ftk,fl)  where  flk,n)  is  given  by  (4.3),  then 


Pp{*CX)  >  rj}  1  1  k  (»!-*!. xi+1)  • 


(4.4) 


n  p4 
i-1  1 


where  g(x1,0, . . . ,0)  -  r^ 
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»  .»  V  ^  /  / 

*  -- 


Vs; 


JV 


Proof .  The  proof  consists  of  finding  a  necessary  and  sufficient 
condition  under  which 


and  then  applying  the  results  of  Pledger  and  Proschan  (1971)  to 

find  the  suprenum  of  the  right  hand  side  of  (4.5)  subject  to 
k  k 

n  p.  -  a.  Clearly,  we  must  have  £  x.  <  n. .  Por  fixed 
i-1  i-1  1 

k  k  k 

l  x.  -  f,  or  equivalently,  for  fixed  ][  (n  -x, )  -  J  n.  ■  f, 

i-1  i-1  1  1  i-1  1 

k  k 

n  (n  -x . )  is  minimized  by  (n.-f)  n  n. .  This  follows  from  the 

i-1  k  i"2 

strict  Schur-concavity  of  II  x,  ,  x  >  0 ,  1  <  i  £  k  (see,  e.g., 

i-1  1  1 

Marshall  and  Olkin  (1979,  p.  78)).  Therefore  we  muBt  have 
x  -  (f , 0 , . . . , 0 ) .  Then  a  necessary  and  sufficient  condition  for 
(4.5)  to  hold  is  that  x  -  (x,0,...,0),  where  x  <  f ’ (k,n)  and 
f 1  (k, n)  is  the  first  positive  integer  f  for  which 


k  k 

(n  -f)  n  n  <  max  II  (n  -x. )  .  (4.6) 

i-2  k  i-1  1  1 

I 

i-1  1 


k 

The  Schur-concavity  of  n  x.  then  gives  (4.3)  as  a  sufficient 

i-1  1 

condition  and  the  subsequent  corollary  gives  a  simple  method  of 
calculating  f'  exactly. 

We  now  assume  that  (4.5)  is  satisfied  and  hence  that 

*  f  k  k  > 

f  (x; a )  -  sup  p  <  l  Y.  >_  n.-x.  +  l  n,»  .  (4.7) 

k  V.i-1  1  A  i-2  ) 

n  p.-a 

i-1  1 


Writing  (4.7)  as  an  iterated  sum  and  noting  that 
I.  (n-x,x+l)  is  a  decreasing  function  of  n  for  fixed  x, 


we  have 


s 


,v, 
V." 
Lv  ; 


where  the  Ui  are  independent  binomial  random  variables  with  para¬ 
meters  Ciij.p^)*  i-2,...,k.  Writing 

Vi  ♦  I  Ui  -  l  I1  Y 
1  i-2  1  i-1  j-1  13 

where  the  are  independent  Bernoulli  random  variables  with 

parameter  p, ,  a  result  of  Pledger  and  Proschan  (1971)  may  be 
1  k  n 

employed  to  show  that  the  upper  tail  of  £  7  Y.  .  is  a  Schur- 

i-1  j-1  13 

convex  function  of  C-lnp1#  -lnp^,...,  -lnp^,  -lnpj* . • . » -lnpj  , . . • 
-lnpk, . . . , -lnpk)  and  therefore  f*(x;a)  -  Ia(n1-x1,  x^+1),  as 

required.  _  _ 

Corollary .  For  each  £,  form  the  vector  z  *  (z^ » *2* *  *  *  * zk3  ^roB 
n  «  (nj.nj, . . . ,nk)  by  continually  reducing  the  maximum  (s)  until 
the  subtractions  total  f+1,  f  >  0.  Denote  by  f  * (k,n^ .nj,  . . .  ,nk) 
the  first  f  for  which 

k  k 

n  z.  >  (n.-f)  n  n.  . 
i-1  1  '  1  i-2  1 

Then  a  necessary  and  sufficient  condition  for  (4.5)  to  hold  is 
that  Xj^  <  f '  (k,n)  . 

Proof.  The  proof  proceeds  exactly  as  for  Theorem  4.1' by  noting 

k  k 

that  2  maximizes  n  r.  subject  to  0  <  r,  <  n.  and  7  r.  « 
k  i-1  1  i  -  i  ^  i 

7  n.-f-l.  This  follows  since  z  is  majorized  by  r  and  the 

i-1  1 

product  is  strictly  Schur-concave . 
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, 


.  .  .  ..............  .w  .  .. .. 


•v. 


(4.3)  reduces  to  (4.1)  which  is 


If  ■  n2  -  . . .  ■  n^, 

Winterbot tom's  (1974)  condition.  However,  s  should  be  replaced 
by  s+1  in  his  formula,  which  also  has  a  sign  error.  As  an 
example,  for  fc  »  2,  n^  »  n2  ■  SO,  from  Winterbottom  (1974),  (4.4) 
is  stated  to  hold  for  x^  <  17  or  nj-Xj  >  S3.  However,  33 >50  <  41 >41 
and  therefore  (4.4)  only  holds  for  Xj  <  13  or  n^Xj  >  37,  as  the 
Corollary  to  Theorem  4.1  shows,  or  the  solution  of  (4.3),  which 
gives  f(2, 50, SO)  -  13.14. 

The  dependence  of  f  on  n  may  be  seen  by  considering  an 
example.  Let  k  ■  2,  n^  ■  S,  n2  «  10.  Then  from  the  Corollary 
following  Theorem  4.1,  (4.4)  only  holds  for  x1  •  0,  whereas  for 
nj  ■  n2  ■  5,  it  holds  for  x^  ■  0,1,2,  and  3.  Thus  the  case  of 
equal  n^,  i«l,2,...,k,  does  not  give  the  minimal  f.  In  fact,  it 
may  be  seen  that  if  n^  >  2nlt  then  (4.4)  holds  only  for  Xj  »  0. 

5.  Concluding  Remarks 

Prom  Table  1,  it  seems  reasonable  to  conjecture  that  (3.5) 
is  valid  for  those  values  of  a,k,n  likely  to  arise  in  practice. 

The  authors  are  continuing  to  investigate  the  problem  and  hope 
to  report  more  precise  conditions  for  the  validity  of  (3.5)  in 
subsequent  work. 
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